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HUBBLE SPACE TELESCOPE SPECTROSCOPY OF BROWN DWARFS DISCOVERED WITH THE 

WIDE-FIELD INFRARED SURVEY EXPLORER 
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Edward L. Wright®, Peter R. Eisenhardt f , M. F. Skrutskie g , Roger L. Griffith®, V Kenneth A. Marsh 1 

ABSTRACT 

We present a sample of brown dwarfs identified with the Wide-fteld Infrared Survey Explorer (WISE) 
for which we have obtained Hubble Space Telescope (HST ) Wide Field Camera 3 (WFC3) near-infrared 
grism spectroscopy. The sample (twenty-two in total) was observed with the G141 grism covering 
1.10—1.70 /tm, while fifteen were also observed with the G102 grism, which covers 0.90—1.10 /am. The 
additional wavelength coverage provided by the G102 grism allows us to 1) search for spectroscopic 
features predicted to emerge at low effective temperatures (e.g. ammonia bands) and 2) construct 
a smooth spectral sequence across the T/Y boundary. We find no evidence of absorption due to 
ammonia in the G102 spectra. Six of these brown dwarfs are new discoveries, three of which are found 
to have spectral types of T8 or T9. The remaining three, WISE J082507.35+280548.5 (Y0.5), WISE 
J120604.38+840110.6 (Y0), and WISE J235402.77+024015.0 (Yl) are the nineteenth, twentieth, and 
twenty-first spectroscopically confirmed Y dwarfs to date. We also present HST grism spectroscopy 
and reevaluate the spectral types of five brown dwarfs for which spectral types have been determined 
previously using other instruments. 

Keywords: stars: low-mass, brown dwarfs: individual (WISEA J032504.52—504403.0, WISEA 
J040443.50—642030.0, WISEA J082507.37+280548.2, WISEA J120604.25+840110.5, 

WISEA J221216.27—693121.6, WISEA J235402.79+024014.1) 


1. INTRODUCTION 

The Wide-fteld Infrared Survey Explorer (WISE) has 
been very successful at identifying the coolest brown 
dwarfs in the Solar neighborhood. WISE provided all 
sky coverage at four mid-infrared wavelengths centered 
at 3.4, 4.6, 12 and 22 /am (Wl, W 2 , IU3, and W 4). 
This coverage has offered the ideal dataset with which to 
identify cool brown dwarfs whose spectral energy dis¬ 
tributions peak at mid-infrared wavelengths. Specifi¬ 
cally, WISE was designed so that the W 1 band coincides 
with a deep water+methane absorption feature and the 
W2 band coincides with a region largely free of opac¬ 
ity. Thus, the WISE W1 — 1+2 color h as been especially 
useful for identifying late type dwarfs (Kirkpatrick et al 


201 1[ Cushing et al. 2011, Kirkpatrick et al.|2012[|M 


lace 


et al.j|2013| and |Thompson et al.||2013|). Seventeen of 

the eighteen spectroscopically coniirmed brown dwarfs 
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WISE ([Cushing 

et al. 2011[ Kirkpatrick et al.|2012, Tin- 

ney et al. 

2012 

, Kirkpatrick et al. 2013, Cushing et al. 

2014a[ and Pinfield et al. 2014). The Y0 dwarf WISE 


dwarf, itself i dentified by WISE (Liu et al. 2012 


et al.| 2014a). Three objects, \VD 0806—661(Luhman 


Leggett 


et al.|l201l|)/CFBDSIR J1458+1013B flLiu et ai.||2011) 
and WISE 0855—0714 ( Luhman||2014 ), likely have ettec 


tive temperatures similar to (or less than) the above Y 
dwarfs, but have yet to be spectroscopically confirmed. 

Ground-based follow-up observations at the mid- 
infrared wavelengths where the spectral energy distribu¬ 
tions of cold brown dwarfs peak are nearly impossible due 
to the high thermal background, forcing follow-up obser¬ 
vations to shorter wavelengths where these dwarfs are ex¬ 
tremely faint. As a result, obtaining moderate signal-to- 
noise (S/N) near-infrared spectra for many of the coldest 
WISE candidates has only been capable with the Wide 
Field Camera 3 (WFC3) aboard the Hubble Space Tele¬ 
scope (HST). HST grism spectroscopy of the latest type 
brown dwarfs has been i nvaluable in the study and classi- 
fication of t hese objects (|Cushing et al.|2011 , Kirkpatrick 


et al. ||2012[ |Kirkpatrick et al.||2(J13[ and |Cushing et al. 
2014a). 


As effective temperatures cool below 600 K, several 
spectroscopic features are predicted by model atmo¬ 
spheres to arise in the Y-band spectral region around 
1.07 /im. These include the emergence of ammonia ab¬ 
sorption components and the disappearance of optical 
alkali resonance lines which have broad wings th at are 
predicted to extend into the nea r-infrared (Burrows et 
al. 2000 Kirkpatrick et al. 2012|). (HST) WEC3 G102 
spectroscopy (0.90—1.10 /ini) allows us to investigate this 
additional wavelength range for such features, as well as 
inspect for differences (and similarities) as a function of 
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spectral type. 

In this paper, we present our HST brown dwarf spec¬ 
troscopic sample, including three new Y-dwarfs: WISEA 
J082507.37+280548.2, WISEA J120604.25+840110.5, 
and WISEA J235402.79+024014.1. We also present im¬ 
proved (higher S/N) HST spectroscopy of five additional 
brown dwarfs for which near-infrared spectroscopy has 
been published previously, and present HST G102 grism 
spectroscopy for a sample of 15 late T and Y dwarfs. We 
first present the six new WISE brown dwarf discoveries, 
then present the new HST WFC3 spectroscopy of previ¬ 
ously identified brown dwarfs. We then present and ana¬ 
lyze the Y-band spectra for the 15 brown dwarfs observed 
with the G102 grism and construct a complete spectral 
sequence across the T/Y boundary. Lastly, we estimate 
physical properties of our entire sample of brown dwarfs 
by atmospheric model fitting. 


2. THE SAMPLE 

The selection criter ia for brown dwarf s in th is study are 
described in detail in Kirkpatrick et al. (2012). A11WISE 
source catalog positions and photometry for all 22 brown 
dwarfs in this study including the six new discoveries are 
given in Table 1. (Note: A detailed analysis of the HST 
spectrum of WISEA J182831.08+265037.6 will be pub¬ 
lished separately in Cushing et al., in preparation). We 
provide data from the A11WISE catalog because of its 
improved astrometric accuracy and increased sensitivity 
in the W1 and W2 bands compared to the WISE All-Sky 
catalog^*] An in depth charact erization of the A1 1 WISE 
source catalog can be found in Kirkpatrick et al. (2014) 
and the A11WISE Explanatory Supplement 10 . Hereafter, 
source names are abbreviated using the first four dig¬ 
its of the right ascension and declination for each object 
(e.g. WISEA J120604.25+840110.5 is WISE 1206+8401). 
Finder charts for each of the new discoveries are pre¬ 
sented in Figures 1—6. 

3. OBSERVATIONS 

3.1. Photometry 

3.1.1. WFC3/HST 

Direct images of each brown dwar f were obtained with 
the Wide Field Camera 3 (WFC3; |Kimble et al.||2008) 
aboard the Hubble Space Telescope. Each grism ob- 
servation (Section 3.2) requires an accompanying direct 
image, necessary for locating sources and determining 
source sizes. Positions and source sizes are then used 
to determine the placement and size of the correspond¬ 
ing grism extraction apertures as well as the wavelength 
zero-point. These images can also be used for photomet¬ 
ric purposes. Direct images for the G141 grism observa¬ 
tions were obtained with either the F140W ( X p = 1392.3 
nm) or F125W filters ( X p = 1248.6 nm), while images 
for the G102 grism observations were obtained with the 
F105W filter ( X p = 1055.2 nm), where X p is the “pivot 
wavelength” (see |Tbkunaga & Vacca||2005). The 1024 x 
1024 HgCdTe detector used by WECJ3 has a plate scale 
of ^0'/13 per pixel, resulting in a total field of view of 
123" x 126". Multiple images were obtained for each 
s ource and were combined using AstroDrizzle (Gonzaga 
eTaLp0l2 ). 


10 See http://wise2.ipac.caltech.edu/docs/release/allwise/expsup/ 



Figure 1. Finder chart for new brown dwarf discovery WISE 
0325-5044. The top panel is the HST WFC3 F125W image, while 
the bottom three panels are WISE channels 1—3, from left to right. 
All panels are centered on the WISE position of the brown dwarf, 
which is also indicated by a red circle. In the top panel, the HST 
brown dwarf position is indicated by a blue circle. North is up and 
East is left in each panel. 
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Figure 2. Same as Fig. 1 for WISE 0404-6420. 

To test the consistency of our measured magnitudes, 
we performed aperture photometry on both the individ¬ 
ual flat field frames (fit - multiplied by the WFC3 IR 
pixel area map to account for geometric distortion) and 
the final drizzled image created by AstroDrizzle (drz). 
In several instances, we found differences between the 
magnitudes derived from the fit and drz ima ges of up 


to ~0.2 mag, similar to the results of [Kir kpatr ic k et al. 
(2012). Further investigation revealed the cause of these 
discrepancies to be liberal default values for the cosmic 
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Figure 3. Same as Fig. 1 for WISE 0825+2805. Figure 5. Same as Fig. 1 for WISE 2212-6931. 
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Figure 4. Same as Fig. 1 for WISE 1206+8401. 


Figure 6. Same as Fig. 1 for WISE 2354+0240. 


ray rejection algorithm of AstroDrizzle. On occasion, 
this algorithm will flag a central pixel of a source as a 
cosmic ray, thus reducing the total flux in the combined 
final drz image compared to the original fit images. Be¬ 
cause cosmic rays are rejected in the calwf3 calibration 
program used to create the fit images, we corrected this 
issue by first turning off the cosmic ray rejection algo¬ 
rithm when combining the fit frames with AstroDrizzle. 
The drz images created by AstroDrizzle in this way were 
found to contain several artifacts (not cosmic rays) that 
are not present in the drz images created by AstroDriz¬ 
zle when the cosmic ray rejection algorithm is turned 
on. We remove these artifac ts by em ploy ing a Lapla- 
cian edge detection algorithm (van Dokkumi2001) to the 


drz images created by AstroDrizzle with the cosmic ray 
rejection algorithm turned off. Aperture photometry is 
then performed on these final images. 

One result of the drizzling process is that the noise in 
adjacent pixels is correlated. For this reason, determin¬ 
ing the rms noise in the local background using common 
background subtraction methods, such as a sky annulus, 
are inapplicable. Instead, we estimate the background 
and its uncertainty by applying the same aperture used 
for the source photometry (O'/d) to 1000 random star-free 
positions (determined via a 3cr clip) on each image. We 
take the mean and standard deviation of these measure¬ 
ments as the background level and its uncertainty. This 
uncertainty is added in quadrature to the signal uncer- 
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tainty, determined by dividing source counts by the gain 
and total exposure time, to determine the final source 
uncertainty. Magnitudes were calculated using the Vega 
system with zeropoints of 25.1845, 25.1439, and 25.4523 
for F140W, F125W, and F105W images, respective!}^ 
HST photometric magnitudes and uncertainties are given 
in Table 2. 


3.1.2. IRAC/Spitzer 

The entire sample was observed wit h the InfraRed Ar¬ 
ray Camera (IRAC; Fazio et al. 2004) aboard the Spitzer 
Space Telescope. Each dwarf was imaged with IRAC 
channels 1 and 2 centered at 3.6 and 4.5 /im (hereafter, 
chi and ch2). A detailed description of the IRAC data 
reduction procedur e used to determine chi an d ch2 mag¬ 
nitudes is given in Kirkpatrick et al. (2011). chi and 
ch2 magnitudes for each dwarf are given m Table 2. 
For WISE 1639—6847, which was blended with another 
source in its chi image, only a ch2 magnitude is given. 


3.2. Spectroscopy 

HST spectroscopic observations were carried out us¬ 
ing the Wide Field Camera 3 (WFC3) as part of Cycle 
18 program 12330 (PI: Kirkpatrick), Cycle 19 program 
12544 (PI: Cushing) and Cycle 20 programs 12970 (PI: 
Cushing) and 13178 (PI: Kirkpatrick) and are summa¬ 
rized in Table 3. The G141 (1.1—1.7 /rm, R = A/AA 
« 130) and G102 (0.9—1.1 /im, R ~ 210) spectroscopic 
re ductions were perform ed using the methods described 
in | Cushing et al.| (|2011 ), with the exception of one final 
step, described in the following paragraphs. 

Because the G141 and G102 grism modes are slitless, 
source spectra are occasionally influenced by photons 
from a nearby source or uneven background fluctuations. 
For this reason, we developed a custom source extrac¬ 
tion routine that allows us to define source apertures 
and background regions on the individual stamp images 
of a spectrum produced with the AXEDRIZZLE routine 
(e.g. Figure 7). Once an aperture is defined, we use the 
aperture corrections of Kuntschner et al. (2011) and the 
appropriate grism sensitivity curve^] to flux calibrate 
the final spectrum. T here is also a 2% abso lute flux 
calibration uncertainty (Kuntschner et al.||2011) which 
we account for in all synthetic photometry calculations 
(Section 5.3) and model fits (Section 5.4). An example 
of the functionality of our extraction technique is shown 
in Figure 7, where the effects of uneven background fluc¬ 
tuations are mitigated. 

Several sources had multiple HST visits (i.e. spectro¬ 
scopic observations occurring on different dates). In prin¬ 
ciple, direct and spectroscopic images from multiple vis¬ 
its can be combined using standard AstroDrizzle routines 
to produce final images with improved S/N. Yet, because 
our sources are nearby brown dwarfs, they tend to have 
significant proper motions, even on time scales as short as 
a month. For HST grism observations, identical positions 
are critical for the placement of extraction apertures and 
wavelength zeropoints. For sources with multiple visits, 
AXEDRIZZLE stamp images from each visit are median 
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Figure 7. Top: Stamp image of G141 the grism spectrum of 
WISE 2056+1459. Solid red lines indicate the defined extraction 
aperture. Dashed lines show the regions used for background fit¬ 
ting/subtraction. Bottom: Comparison of the AXEDRIZZLE ex¬ 
tracted spectrum of WISE 2056+1459 (black) and the spectrum 
extracted with our routine (red). 


combined to produce a final spectroscopic image. Spec¬ 
tra are then extracted using the method described above. 
Objects with both G141 and G102 data are stitched to¬ 
gether at 1.1 /im with no scaling done between the two 
spectra. The complete spectroscopic sample is shown in 
Figures 8 and 9. We note here that the dep ressed H 


band peak of WISE 0647—6232 seen in Kirk patrick et al. 
(20l3|) does not appear in our reduced spectrum. This 
is likely due to either the added signal of our spectrum 
or our improved spectral extraction technique. WISE 
0535—7500 is located in an extremely crowded field, and 
thus retrieving a clean G102 spectrum was particularly 
difficult, even with multiple observations at multiple roll 
angles. While we did manage to extract a G102 spectrum 
for this object, we urge caution in its interpretation, as it 
is still likely contaminated. We also note that a different 
roll angle allowed us t o extract a cleane r J-ba nd spec¬ 


trum than that found in Kirkpatrick et al. (2012), though 


the zeroth order light from a nearby object fell directly 
in the H-band portion of WISE 0535—7500’s spectrum. 
For this reason, we do not use the H-band portion of this 
object’s spectrum. 

Upon inspection of their spectra, two objects that met 
the selection criteria were found to not be brown dwarfs. 
The spectra for both of these objects were flat and devoid 
of spectroscopic features and are likely extragalactic in 
nature. The A11WISE positions and photometry for these 
contaminants is provided in Table 1. 


4. ANALYSIS 


11 http://www.stsci.edu/hst/wfc3/phot_zp_lbn 

12 http://www.stsci.edu/hst/wfc3/analysis/grism_obs/wfc3- 
grism-resources.html 


4.1. New Discoveries 

Spectral types for the six new brown dwarf dis¬ 
coveries (WISE 0325-5044, WISE 0404-6420, WISE 
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Figure 8. All HST spectra (black) ordered by spectral type, normalized at 1.27 fi m, offset by constants, and compared with corresponding 
spectral standards (red). The spectral standard s and the instrume nts from which the comparison spectra were ob tained are as follows : 
T8-2MASSI J0415195—093506 (IRTF/SpeX; Burgasser et al. 2006), T9-UGPS J072227.51-054031.2 (IRTF/SpeX;|Cushing et al. 2 0TTt, 
Y0—W ISE J173835.53+273259.0 (HST/WFC3; Cushing et al. 2011), and Yl-WISE J035000.32-565830.2 (HST/WFUTf KirkpatricFetl 
|al.|2012| ). 



Wavelength (/tm) Wavelength (urn) 


Figure 9. All HST spectra (black) ordered by spectral type, normalized at 1.27 //in, offset by constants, and compared with corresponding 
spectral standards (red). The spectral standards and the instruments from which the comparison spectra were obtai ned are as follows: 
T8-2MASSI J0415195—093506 (IRTF/SpeX; IB urgasser et al.12 006), T9-UGPS J072227.51-054031.2 (IRTF/SpeX: ICushi ng et al.l20lU, 
Y0—W ISE J 173835.53+273259.0 (HST/WFC3; Cushing et al. 2011), and Yl-WISE J035000.32-565830.2 (HST/WFC3; Kirkpatrick et 
|al.|2012) . 
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0825+2805, WISE 1206+8401, WISE 2212-6931, and 
WISE 2354+0240) were de ter mined following the meth - 
ods of|Cushing et al.|(|2011| and |Kirkpat rick et al.|(|2012|). 
These discoveries are classified via by-eye comparisons of 
the width of the J-band peak to the spectral standard s 


( 2011 ) 


defin ed in|Burgasser et alT| (|2 006|) , |Cushing et al. 
and Kirkpatrick et al. (2013|. Comparisons are snown m 
Figures 8 and 9, while the spectral types are provided in 
Table 1. 

The spectrum of WISE 0325—5044 is an excel¬ 
lent match to t he T8 spectral stan dard 2MASSI 
J0415195—093506 (|Burgasser et al.|2006|), hence we clas¬ 
sify it as a T8. Similarly, WISE 0404—6420 and WISE 
2212—6931 show excellent agreement with the T9 spec¬ 
tra l stan dard UGPS J072227.51-054031.2 (|Cushing et 
al.||2011) and are classified as such. WISE 1206+8401 
is an almost perfect match t o the Y0 spectral st andard 
WISE J173835.53+273259.0 ([Cushing et al.|20TT|, there¬ 
fore we are confident in assigning it a Y0 spectral classi¬ 
fication. Even though the spectrum of WISE 2354+0240 
is somewhat noisy, it shows a J-band peak significantly 
narrower than the Y0 spectral standard and matches 
well with the Y1 standard WISE J035000.32—565830.2 


(Kirkpatrick et al. 2012). Therefore, we assign a spec¬ 
tral type of Y1 for this object, making it the fourth 
such object to receive such a designation, along with 


WISE J035000.32—565830.2 (Kirkpatrick et al 


WISE J064723.23—623235.5 ([Kirkpatrick et al 
and WISE J154151.65-22502 


2013), 

9 (see Section 4.2 + The 


2012 ), 


most ambiguous of all classifications was that of WISE 
0825+2805, which is shown separately in Figure 10. A 
comparison of WISE 0825+2805 with the Y0 and Y1 
spectral standards shows that the width of the J-band 
peak for this object is intermediate between the two. For 
this reason we assign WISE 0825+2805 a spectral type 
of Y0.5. 

While all of these discoveries are part of an extensive 
parallax program to measure their distances, we provide 
preliminary estimates here based s olely on W2 photom- 
etry and the polynomi a l fits from Dupuy & Liu (2012) 
and Kirkpatrick et al. (2012). Distances estimates are 
provide in Table 4. Uncertainties are based solely on the 
photometric uncertainty of the W2 measurement. 

4.2. Improved Spectroscopy of Previously Identified 
Brown Dwarfs 

Figures 8 and 9 also show HST spectroscopy of 
five previously known brown dwarfs, including WISE 
J033 5+4310 (T9-|Mace et al.|201 3D, WISE J07344-7157 
(YO- pGrkpatriclt 


et al.l [20l2]r WISE J1541-2250 


(“Y0.57” - Kirkpatrick et al. 2012), WISE J1639-6847 
(Y0-Y0.5 -jTinney et al.n2UT2TTmd WISE J2220-3628 


(Y0- Kirkpatrick et al. 2012). Spectral types were re¬ 
examined using the methods described in Section 4.1. 
Spectral types determined in t his way agree well w i th the 


spectral types determined m_ Mace et al i ^_ 

Patrick et al. (2012), and Tinney et al. "*( 2012) 


(2013]) 

2012 ). 


Kirk- 
WitlT 

the improved "S/JN from HSl\ we reclassify W1541—2250 
as Y1 because the width of its J-band peak is a much 
better match to the Y1 standard than that of the Y0 
standard. This determination super sedes the “Y0.5?” 
suggested in Kirkpatrick et al. (2012). The J-band peak 
of WISE 1639—6847 matches well with the Y0 standard, 
in agreement with the Y0-Y0.5 spectral type determined 



Figure 10. Spectral classification of WISE 0825+2805. WISE 
0825+2805 is shown in black, while spectral standards are shown 
in red. 


m 


Tinney et al. (2012). However, the Y-band peak and 


Y — J color of this object is unusual compared to other 
Y0 dwarfs (see Section 5.1 and Figure 13), therefore we 
classify WISE 1639—6847 as YOpec. 


5. DISCUSSION 
5.1. Spectral Sequence 

In Figure 11 we show the entire combined HST 
G102/G141 near infrared spectra for a complete sequence 
of late-type brown dwarfs from T8 to Yl. There is a 
distinct trend for the Y, J, and H band peaks as temper¬ 
atures decrease. The Y band peak appears more sym¬ 
metrical going from T8 to Yl. This occurs because the 
blue wing of the Y-band peak becomes more suppressed 
for later spectral types. The wavelength at which the 
peak reaches a maximum also occurs at a different loca¬ 
tion as temperature decreases, gradually becoming bluer 
with spectral type. The J-band peak becomes narrower 
as temperatures decrease, while the H-band peak appears 
to become more symmetrical for later spectral types. The 
peaks of both the J- and H-bands occur at the same po¬ 
sition throughout the sequence. 

G102 observations were made for five of the six 
new brown dwarfs presented in Section 4.1 (the ex¬ 
ception being WISE 2354+0240), all five dwarfs pre- 



WISE 1639—6847 shows an increase in the relative 
height of the Y-band peak compared to other Y0 type 
dwarfs. The peak is also shifted blueward from the peak 
of the T9 standard UGPS J072227.51-054031.2. This is 
a surprising result considering that we find the J-band 
peak matches well with the Y0 standard (see Figure 9). 
Because the Y-band shape of this object is so unusual, 
we classify it as a YOpec. We note that a similarly large 
Y-band peak is seen in the FIRE spectrum of WISE 













































































































HST Spectroscopy of WISE Brown Dwarfs 


7 



Figure 11. The complete spectral sequence from T8 to Yl. 


1639—6847 presented in Tinney et al.j (2012) (see their 
Figure 7). 

The YOpec dwarf WISE 1405+5534 (which is catego¬ 
rized as ‘pec’ due to a wavelength shift of its H-band 
peak) shows a Y-band peak slightly different than that 
of the T9 and Y0 dwarfs. Though noisy, the Y-band peak 
of WISE 1405+5534 is more symmetrical about its max¬ 
imum than the normal T9-Y0 dwarfs. We also note here 
that the J-band peak appears slightly narrower than that 
of the Y0 standard (see Figure 9). For this reason, we 
reclassify WISE 1405+5534 as Y0.5. The near-infrared 
colors of WISE 1405+5534 support this change (see Sec¬ 
tion 5.3). 

WISE 2209+2711 (Y 0 :-|Cushing et al. 2014a|), while 
noisy, shows a very peaked appearance, more similar to 
that of the Yls than the YOs. This may be providing 
a hint that this object is possibly Yl or intermediate 
between Y0 and Yl. 

5.2. Nm 


The appearance of ammonia (NH 3 ) in the near infrared 
spectra of cool brown dwarfs has been suggested for some 
time as the harbinger of the Y spectral class ([Burrows 


et al. 
Mode! 


2003 Leggett et al. 
s predict that 1 NH 3 a 


2007 and Kirkpatrick 2008p 
Isorption bands could become 


prominent in brown dwarfs with temperatures <450 K. 


Cushing et al. (2011) show a possible detection of NH 3 on 
the blue wmg of the H-band peak of WISE 1738+2732, 
but this NH 3 feature is itself confused with a strong H 2 O 
band. An additional NH 3 absorption component is pre¬ 
dicted to emerge that will affect the shape of the Y-band 
spectrum of extremely cool dwarfs, eventually causing 
a bifurcation of the peak around 300 K. So far, a firm 
detec tion of NH 3 absorpt io n at 1.03 im 1 has bee n elu¬ 
sive (Leggett et al. 2013). Leggett et al. (2013) have 
suggested that vertical mixing can explain the lack of 
NH 3 in the spectrum of the Y0 dwarf WISE 2056+1459. 
Vertical mixing keeps the atmospheric chemistry from 
reaching equilibrium, thus enhancing the N 2 ab undance 

and decreasing the NH 3 abundance (Saumon et al. 2003 

1 1 / R — '+1 


Hnbeny fe Burrows|20d7 , Morley et al.|2014| ). Mow, with 
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Figure 12. HST G102 spectra for WISE brown dwarfs (black). Each spectrum is normalized at 1.27 n m (the J-band peak) and ordered 
sequentially by spectral type from top to bottom, then left to right. For comparison purposes, the spectrum of the T9 standard UGPS 
J072227.51—054031.2 is included in red on each plot. 


an increased sample size of late-type brown dwarfs for 
which Y-band spectra are available, we can investigate 
the effects of NH 3 in these atmospheres. If the lack of an 
NH 3 detection in the spectrum of the Y0 dwarf WISE 
2056+1459 is due to vertical mixing, as suggested by 


Leggett et al. (2013), then such a phenomenon may be 
ubiquitous in early-Y dwarfs. We find no obvious evi¬ 
dence of the presence of NH 3 in the Y-band spectra for 
any of our late-type dwarfs. 


5.3. Colors 

We synthesized Y —, J —, and H —band photometry us¬ 
ing the MKO-NIR photometric system for each brown 
dwarf in our sample. We also synthesized HST F105W, 
F125W, and F140W photometry using available system 
throughput tableful Uncertainties are computed using 
a Monte Carlo approach. For each spectrum, a value is 


selected from a normal distribution using the flux den¬ 
sity at each wavelength as the mean and the flux density 
uncertainty as the standard deviation. Synthetic pho¬ 
tometry is performed on the resulting spectrum, and the 
process is repeated 1000 times. We take the mean pho¬ 
tometric value as the true value, and the standard de¬ 
viation of the distribution as the uncertainty. The re¬ 
sulting magnitudes and uncertainties are listed in Ta¬ 
ble 5. To investigate trends as a function of spectral 
type, we also synthesize photometry on the entire sam¬ 
ple of T dwarfs available from the Spex Prism Spectral 
Librar}p^[as of June, 2014) using the same method de¬ 
scribed above. Colors as a function of spectral type are 
shown in Figure 13. 

Our synthetic YJH photometry generally agre es we ll 
with published ph o tomet r y from Cushing et al . ( 2011 |), 
Kirkpatrick et al. (2012), IVlorley et al. Q 2012 [ ) , bupuy 


13 http://www.stsci.edu/hst/wfc3/ins_performance/throughputs 


14 http://pono.ucsd.edu/~adam/browndwarfs/spexprism/ 
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& Krausl (|2013|), [Leggett et al. 
(|2014MC ushing et al. (201^7 


(2013| 

ilrkpal 


Beichman et al. 


irkpatrick et al. 


(2014) 

ons(51 


and|Leggett et al.|(|2(JMb|), with only a tew exceptions (56 
out of 62 (87%) measurements have differences that are 
within 3cr of zero). Our measured H magnitude for WISE 
0410+1502 and J magnitude for WISE 1405+5534 dif¬ 
fer signifi cantl y from the measured valu es in |Cushing et 


al. (2011) and Kirkpatrick et ah J 2012|, though th ey do 


agree with the values quoted in Leggett et al. (2013). Our 
J- and H-band magnitudes for WI SE 1541—2 250 disagree 


with the valu es given in Morley et al. (2012) and Leggett 


et al. (2014b), respectively. U sing our measured magni¬ 


tudes, WISE 1541—2250 does not stand out prominently 
from the other Y1 dwarfs in color-color and color ver¬ 
sus spectral type diagrams. Our measured J magnitude 
for WISE 1738+2732 disagrees with the value quoted in 


Leggett et al. ( 2013), though does a gree with the mea¬ 
surements from | Cushing 
ly, our m< 


et al. 


(2011) and Kirkpatrick et 


ah] p012j) . Lastly, our measured J and H magnitudes lor 
WIS E 2056+1459 agree w ell with the measured value s 


IM2I) 

,t et al.l 


from Cushing et al. (2011) and Kirkpatrick et al. 

though do not agree with the values from |Leggett _ 

(2013). While we cannot rule out extreme variability 
as the source of magnitude differences, it is more likely 
that some measurements are erroneous. Since our pho¬ 
tometry comes from our slitless grism spectroscopy, one 
potential source of error in our measurements is contami¬ 
nation from nearby bright objects or uneven background 
fluctuations, though our reduction process (Section 3.2) 
should minimize such effects. 

Previous investigations of color/spectral typ e relat ions 

y (Lodieu"etl 

Leggett et al.||2013[ iLeggett et al.||2(J14bh . We 
see similar changes in Figure 13, namely, the Y — J (and 
F105W — F125W) color plummets for spectral types 
later than T8, and the J — H color plateaus at ^T7 
and slowly turns redder for the Y dwarfs. As noted by 
Leggett et al. (2010|, the blueward trend in the Y — J 
colors of late type dwarfs may be due to K I condens¬ 
ing into KC1 as temperatures decrease. This condensing 
weakens the broad 0.77 /im K I doublet, resulting in a 
brighter Y-band. A similar tre nd i s seen the the z — J 
colors of Y dwarfs in Lodieu et al. (2013). The redward 


al 


trend of the J — H color Tor the latest Y dwarfs occurs 
because the J-band peak for these objects continues to 
narrow, while the H-band peak does not (see Figure 11). 

Figure 14 shows the Y — J (when available) versus 
J — H and J — H versus J — W2 colors for our complete 
sample of brown dwarfs along with synthesized colors 


from the mo del s pectra of Saumon_et al. (2012), Morley 

1 


et al 


(20121), and 


Morley et al. J 2014[ ). While neither the 


cloudy or cloud-free models reproduce the trends seen for 
our measured dwarfs, the cloudy models do a better job 
approximating the colors of our late-type d warf s ample 
than the cloud-free models. Morley et al. (2012) show 
that the inclusion of clouds is an efficient way to redden 
the J — H colors of late type brown dwarfs. Figure 14 
also shows that models with lower surface gravities and 
lower f se d values do not fit the observed color trends as 
well. 


5.4. Physical Parameters 

The effective temperatures of brown dwarfs are typi¬ 
cally estimated in two ways. If the distance to a brown 


dwarf is known and a significant fraction of its emergent 
flux can be measured, then its bolometric luminosity can 
be computed. Since brown dwarfs all have similar radii 
due to the equation of state of their partial degenerate 
interiors, an effective temperature can then be c omp uted 


using the Stefan Boltzmann law (e.g. Golimowski et al. 
2004, Vrba et ah||2004| Dupuy fc Kraus||2013| ). Alter¬ 


natively, model spectra can be compared directly to ob 
served spectra and/or photometry to infer not only the 
effective temperatures of brown d warfs, b ut also their 
surface gravity an d metal licity (e.g. Cushing et al.||2008 


Witte et al. 2011, and Leggett et al. "2013K Because 
accurate distances to many of the dwarfs in our sam¬ 
ple are not known, we have taken the second approach 
and compared our near-infrared spectra and mid-infrared 
photometry to atmospheric models in order to estimate 
their atmospheric parameters. 

In particular, we compare are spectra and photometry 
to several different sets of solar metallicity modeEj in¬ 
cluding the cloud-fre e models from Saumon et al. (2012), 
the cloudy models of |M orley et al. 2012]) , and the cloudy 
models of|lVlorley et al.| (|2014|). The latter set is a new 
suite of models that in addition to including sulfide and 
chloride condensates, also includes water ice clouds for 


obj ects w ith T e ff < 450 K. Some models from Morley et 


al. (2014) are also “partly cloudy”, calculated with 50% 


Houdy" and 50% hole surface coverage. The models are 
arranged in a grid of T e //, log g, and f se d values. The 
f se d parameter describes the efficiency of cloud sedimen¬ 
tation in the models, where larger f se d values correspond 
to larger grain sizes, which lead to an increased sedimen¬ 
tation efficiency, and hence thinner clouds. W e adop t 


the fsed = “nc” or “no clouds” notation of Cushin g et 


al. (2008) to refer to the cloud-free models. Table 6 sum¬ 


marizes the T eff, log g , and f se d combinations that were 
used for model fitting. 

In order to compare the models to our observed spec¬ 
tra, each model spectrum was resampled to be uniform in 
ln(A) space, smoothed with a Gaussian kernel to a resolv¬ 
ing power of HST WFC3, and resampled onto the same 
wavelength grid as the HST WFC3 spectrum. Synthetic 
flux density values in the [3.6] and [4.5] Spitzer bands 
were computed from th e models as described in Cushing 
et al. (2006). Following Cushing et al. ( |2008 ), we evalu¬ 
ate a goodness of fit parameter for each model spectrum, 
defined as 


G k = Y 


i=1 


fi CkJ~k,i 


(1) 


where n corresponds to each data pixel, fi is the flux 
density of the data, Fk,i is the flux density of the model 
and cq are the errors for each observed flux density. 
The Ck parameter is a multiplicative constant given by 


= E/VmM 

E-TmM 2 ' 


( 2 ) 


Ck is equivalent to (R 2 /d 2 ), where R is the stellar radius 
and d is the distance. For each model we compute a Ck 
value and corresponding Gk value and identify the best 
fitting model as that model with the lowest overall Gk 
value. Because in most instances the difference between 
the best fitting models for each dwarf are marginal, we 
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Figure 13. Left: Y — J color as a function of spectral type for T and Y dwarfs. T dwarfs from the SpeX Prism Spectral Library are 
represented by blue triangles, while synthesized colors of brown dwarfs in this study are represented by cyan circles. Small offsets have 
been added along the abscissa for differentiation purposes. Center: J — H color as a function of spectral type for T and Y dwarfs. Right: 
HST colors as a function of spectral type for T and Y dwarfs. Synthetic photometry of brown dwarfs in this study are represented by cyan 
circles, while colors found using aperture photometry are represented by orange squares. 


proved a range of best fit model parameters from the 
three best fitting models for each brown dwarf in Table 
7. The ( R 2 /d 2 ) and y 2 values in Table 7 refer to the 
model with the best fit. 

The derived effective temperatures are generally cor¬ 
related with spectral type in that the effective temper¬ 
atures decrease with increasing spectral type: WISE 
0325—5044 (the lone T8) has an effective temperature 
range of 550—600 K, the T9 dwarfs have temperatures 
between 500 and 600 K, the T9.5 dwarfs have tempera¬ 
tures between 450 and 500 K, the Y0 dwarfs have tem¬ 
peratures between 400 and 450 K, the Y1 dwarfs have 
temperatures between 300 and 400 K. Exceptions include 
the Y0: dwarf WISE 2209+2711 (T e// = 500-550 K) 
and the Y1 dwarf WISE 0535—7500 which is located in 
an extremely crowded field and thus has contamination 
issues (see Section 5.1). 

A comparison with the effecti ve temperatures derived 
from bolometric luminosities in Dupuy & Kraus (2013) 
find some agreement, though there are several discrep- 
ancies. The main discrepancies app ear for t he ob¬ 
j ects w ith very uncertain parallaxes in|Dupuy & Kraus| 
(WISE 0359-5401, WISE 0535-7500, WISE 
L—2250, and WISE 1639—6847). For these objects 
with uncertain parallaxes, two (WISE 0359—5401 and 
WISE 16 39—6847) are found to have T e ff values much 
colder in Dupuy & Kraus (2013) than our model-fitting 
derived values, while the other two (WISE 0535—7500 
and WISE 1541—2250) are found to be much warmer. 
It is unclear whether these discrepancies are solely due 
to parallax uncertainties or other systematic errors (see 
below). 

Although atmospheric models can be computed for any 
given T e ///log g values, not all combinations are phys¬ 
ical because theoretical evolutionary models limit the 
T e jj/log g values that brown dwarfs can obtain. We 
therefore compared our derived T e /j/log g values to the 
low- mass, solar metallicity, clou d-free evolutionary mod¬ 
els of Saumon & Marley (2008) in order to confirm that 
our best fit models are all physical (Figure 15). We note 
that the atmospheric models used in our analysis were 
not used as boundary conditions for the evolutionary 
models and thus our comparisons are not technically self 
consistent. Inspection of Figure 15 suggests that several 


is 


of the derived surface gravities are not physical (those 
objects with log g values of 5 and T e ff estimates < 400 
K). Evolutionary models can also be used to estimate 
the ages of brown dwarfs. Most of our sources have ages 
of several Gyr, as expected for objects near to the Sun, 
though some objects were found to have younger (a few 
hundred Myr) or older (> 5 Gyr) ages than expected. 
Indeed, a trend is apparent whereby the earlier spectral 
types in our sample (T8 and T9) are preferentially fit 
with lower log g values and hence younger ages, while 
the latest spectral types (Y0.5 and Yl) are preferentially 
fit with larger log g values (some of which are unphysi¬ 
cal). Although it is possible this trend is real, it is more 
likely a result of systematic errors in the model fits. Fits 
with unphysical gravities and are indicated as such in the 
age column of Table 7. 

While the derived temperatures are generally consis¬ 
tent with previous work, the fits to the data are, in many 
cases, rather poor. The discrepancies between the data 
and models generally fall into two categories. Because 
the HST spectra contain over two hundred data points, 
the two IRAC points carry little weight in the fitting 
process and thus the models “match” the near-infrared 
data reasonably well but often miss the Spitzer /IRAC 
points badly. Figure 16 shows the best fits to the T8 to 
Yl spectral sequence from Figure 12. In some instances 
(e.g. WISE 0325-5044, WISE 1541-2250), the Spitzer 
model fluxes can be off by factors of 0.2 to 4, respec¬ 
tively. Since much of the energy of cold brown dwarfs 
emerges at 5 /im, the bolometric luminosity of the best¬ 
fitting models disagree with the observations by similar 
factors. A factor of 2 to 5 in L^i can result an effective 
temperature difference of up to 200 K (see Figure 15). 

Closer inspection also reveals large differences between 
the models and observations in the shapes of the near- 
infrared bands. In particular, the heights of the model 
i7-band peaks are significantly different than the ob¬ 
servations, with the observed H-band peak consistently 
being higher then that of the best-fitting model. The 
height, width, and position of the J-band peak is gen¬ 
erally well matched to the models. For the two latest 
spectral types, the J-band peak of the observations does 
not reach the same heights as those from the best-fitting 
models. The Y-band peaks for the two T dwarfs match 
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Figure 14. Left: J — H versus Y — J s ynth esized colors for brow n dwarfs from our sample, cloud-free |Saumon et al. (2012]) models, and 
clou dy models from Morley et al. (2012) and Morley et al. (2014). Right: J — W 2 versus J — H for brown dwarts trom our sample with 
the|Saumon et ah|([2012 ), Morley et al. (2012), and Morley et al. (2014| models. f se g values apply to the cloudy and partly cloudy models 
( not the cloud-tree models) m all panels. 


well with the models, while the fits to the Y-band peaks 
of the Y dwarfs do not. The main discrepancy in the 
Y-band spectra of the Y dwarfs is the absence in the 
observations of the NH 3 absorption feature expected to 
occur for temperatures <500 K. The problems with these 
particular model fits are nicely summarized by the color- 
color diagrams shown in Figure 14. 

Some of the issues in the near-infrared are likely a re¬ 
sult of the fact that the models do not account for non- 
equilibrium NH 3 chemistry due to vertical mixing in the 
atmosphere (see Section 5.2). Figure 17 shows the near- 
infrared spectrum of WISE 1541—2250 (Yl) along with 
the cross section spectrum of NH 3 . Ammonia absorption 
shapes the entire near-infrared spectrum - from the blue 
wing of the 77-band peak, to the width of the J-band 
peak, to the shape of the Y-band peak. Morley et al. 
(2014) show that when all else is equal, the decrease in 
the NH 3 abundance broadens the J- and 77-band peaks, 


and removes the strong NH 3 absorption feature in the 
Y band. The inclusion of non-equilibrium chemistry, 
which would result in significant changes to the model 
spectra would no doubt result in different effective tem¬ 
perature and surface gravity estimates. Based on the 77 
band alone, the change would be 300 K at a fixed grav¬ 
ity. Given the significant mismatches identified between 
these model spectra and data, we conclude that we can¬ 
not reliably derive reliable atmospheric parameters by 
fitting these model spectra to the observed spectra and 
photo metry of o ur samp le. 


Leggett et al. (2014b) state that the Y, 77, and chi 
model fluxes are too low by about a factor of two com¬ 
pared to photometric measurements for late-type brown 
dwarfs. By extension, this suggests that the J and W2 
magnitudes are not affected and thus the J — W2 color 
may be able to provide a reasonable effective temper¬ 
ature estimate when compared with colors from model 
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Figur e 15. The best fitting model parameters from Table 7 compared with the cloud-free evolutionary models from |Saumon Sz Marleyj 
(|2008|). Lines of constant mass are plotted in solid black for 0.005, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06 Mq from bottom to top. Lines of 
constant luminosity are plotted in solid red for -6.4, -6.0, -5.4, and -5.0 log(L /Lq) from left to right. Isochrones are plotted as dotted blue 
lines for 10, 4, 2, 1, 0.4, 0.2, 0.1, 0.04, 0.02, and 0.01 Gyr from left to right. Small offsets have been added along the abscissa and ordinate 
for differentiation purposes. 


spectra. Figure 18 shows the J — W2 color as a func¬ 
tion of effective temperat ure for three dif f erent models 
- a cloud-free model from Saumon et al. (| 2012 |) (log q 


= 5), a 50% cloud coverage model from IVlorley et al. 
( 2014| ) (log g = 5, f se d = 5, with iT 2 Q ), and a 100 % 


cloud coverage model from|Morley et al.| ( 2012 ) (log g = 
5, fsed = 3, without H 2 O). Because our grid of models 
is, as of yet, incomplete (see Table 6 ), we show the only 
50% cloud coverage model (with H^O ) for which the en¬ 
tire 300—1000 K temperature range is available and the 
100 % cloud coverage model does not extend to tempera¬ 
tures < 400 K. Even when limited to these three models, 
Figure 18 shows that effective temperature estimates for 
J — W 2 values less than typically have ranges greater 
than 100 K. For J — W2 values greater than ~7, varying 
percentages of cloud cover result in minimal differences 
in effective temperature, though the 100 % cloudy model 
does not extend to this range. The range of tempera¬ 
tures for a single J — W2 value will only increase with 
the inclusion of additional models with varying log g and 
f sed values. 

Furthermore, as noted previously, the J—band flux will 
also be affected by the reduced amount of NH 3 present in 


the atmosphere du e to vertical mixing (not only Y and 
77). Figure 12 of Morley et al.| ((20146 displays equilib¬ 
rium and disequilibrium models at 450, 300, and 200 K. 
For the 450 K models, including disequilibrium chemistry 
will change the Y —, J —, and 77— band magnitudes by 
^ 0 . 2 , ^0.3, and ~ 0 . 6 , respectively, implying that these 
near-infrared bands are affected by similar amounts when 
vertical mixing is included. Because of the large tem¬ 
perature ranges for individual J — W2 values, uncertain 
temperature ranges for J — W2 colors due to of our in¬ 
complete model grid, and the similar influence of verti¬ 
cal mixing on the J—band magnitude when compared 
to the Y— and 77— bands, we conclude that estimating 
temperatures based solely on model J — W2 colors is un¬ 
sound. Since our observations are limited to near- and 
mid-infrared wavelengths and the current models do not 
include non-equilibrium chemistry, we lack a method to 
determine accurate effective temperatures for the objects 
in our sample. 

Given the limitations of these particular models, de¬ 
termining effective temperatures by computing bolomety 
ric luminosities may therefore be preferable (e.g. Dupuy 
& Kraus 2013). However, determining luminosities m 
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Wavelength {pm) 


Figure 16. The T8 to Y1 spectral sequence (black) with best fit models (red). The best fitting model parameters (T e ff in K, log g in cm 
s — f se di and the percentage of cloud cover) are given under each object name. The black symbols represent the Spitzer IRAC photometry 
from Table 2, while the red circles indicate the synthetic photometry of the model in the chi and ch2 bandpasses. Error bars for the IRAC 
photometry are smaller than the symbol size. 


this way currently relies on models to estimate the_flux 


levels for up to 50% of the total lum inosity (Dupuy & 
Kraus||20l3 ). Dupuy & Kraus (2013) note that all cur- 


rent atmospheric models provide similar bolometric cor¬ 
rections in the wavelength ranges where flux has not yet 
been directly measured. Nevertheless, minimizing the 
dependence of the derived bolometric luminosities on the 
models should be a high priority, and therefore determin¬ 
ing accurate effective temperatures for the lowest mass 
brown dwarfs should focus on obtaining flux measure¬ 
ments at wavelengths longer than ^5 /im. The near- 
infrared (NIRSpec and NIRCam) instruments and in par¬ 
ticular the mid-infrared instrument (MIRI) aboard the 
forthcoming James Webb Space Telescope could prove in¬ 
valuable in providing the additional flux coverage needed 
to determine accurate effective temperatures for the low¬ 
est mass brown dwarfs. 


6. CONCLUSIONS 

We have presented HST near-infrared spectra for a 
sample of twenty-two brown dwarfs, six of which are new 
discoveries. Three of the new discoveries are classified 
as having spectral type Y, bringing the total number 
of spectroscopically confirmed Y dwarfs to twenty-one. 
Theoretical spectra for the lowest temperature brown 
dwarfs do not yet accurately reproduce the spectra from 
these observations. While disagreements between mod¬ 
els and observations have been known to exist for these 
objects, the inclusion of HST G102 spectra has made dis¬ 
crepancies even more apparent. The primary disagree¬ 
ments occur in regions where models predict^ d etecta ble 
NHs. These spectra, as wel l as those in ILeggett et al. 


(2013), Leggett et al. (2014a), and|Leggett et al. (|2014bf 
show that the strong NH 3 absorption features expected 
to occur in the Y- and H-band spectra of Y dwarfs are 
not present. While the non-equilibrium chemistry for T 
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Figure 17. The near-infrared spectrum of WISE 1541—2250. The 
cross section spectrum of NH 3 at T = 600 K and P = 1 bar (Richard 
Freedman, priv. comm.), indicating the location of prominent NH 3 
absorption bands, is shown in blue. The NH 3 cross section spec¬ 
trum is normalized at 1.033 fi m, inverted, and offset for clarity. 



Figure 18. The J — W 2 color as a fu nctio n of effective tem pera- 


ture for a cloud-free model model from [Saumon e t al. 


= 5), a 50% cloud coverage model from A lor ley et al. 


2012 ) (log g 
(2014) (log g 


= 5, = 5 , with H 9 O), and a 100 % cloud coverage model from 

|Morley et al. (2012) (log g = 5, f se d = 3, without H 2 O). Each 
horizontal gray line indicates a single J — W 2 value for a brown 
dwarf in our sample. 


dwarfs rep roduce observations quite well (Hubeny & Bur- 
rows|2007 ), the latest non-eq uilibrium models do not yet 
agree with the observations (Morley et al. 2014|). That 
the models diverge from accurately describing the spec¬ 
tra of objects at such low temperatures begs for new 
models that include improved non-equilibrium chemistry, 
or the identification of another method for reducing NH 3 
abundances. 
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Table 1 

A11WISE Photometry 


A11WISE Name 

Spec. 

W1 

W2 

W3 

W1-W2 


Type 

(mag) 

(mag) 

(mag) 

(mag) 


New Discoveries 


WISEA J032504.52—504403.0 

T8 

18.430 ± 0.258 

16.209 ± 0.145 

>12.918 

2.22 ± 0.30 

WISEA J040443.50—642030.0 

T9 

18.442 ± 0.178 

15.726 ± 0.063 

>12.977 

2.72 ± 0.19 

WISEA J082507.37+280548.2 

Y0.5 

>18.444 

14.578 ± 0.060 

>11.660 

>3.87 

WISEA J120604.25+840110.5 

Y0 

>18.734 

15.058 ± 0.054 

>12.536 

>3.68 

WISEA J221216.27—693121.6 

T9 

17.259 ± 0.122 

14.873 ± 0.061 

>12.621 

2.39 ± 0.14 

WISEA J235402.79+024014.1 

Y1 

>18.263 

15.007 ± 0.085 

>12.278 

>3.26 


Other Brown Dwarfs 

in this Study 



WISEA J033515.07+431044.7 

T9 

>18.652 

14.515 

± 

0.055 

>11.901 

>4.14 

WISEA J035000.31—565830.5 

Y1 

>18.699 

14.745 

± 

0.044 

12.325 ± 0.282 

>3.95 

WISEA J035934.07—540154.8 

Y0 

>19.031 

15.384 

± 

0.054 

>12.877 

>3.65 

WISEA J041022.75+150247.9 

Y0 

>18.170 

14.113 

± 

0.047 

12.314 ± 0.500 

>4.06 

WISEA J053516.87—750024.6 

>Y1 

17.940 ± 0.143 

14.904 

± 

0.047 

>12.349 

3.04 ± 0.15 

WISEA J064723.24—623235.4 

Y1 

>19.539 

15.224 

± 

0.051 

>12.961 

>4.32 

WISEA J073444.03—715743.8 

Y0 

18.749 ± 0.281 

15.189 

± 

0.050 

>12.959 

3.56 ± 0.29 

WISEA J094306.00+360723.3 

T9.5 

18.176 ± 0.297 

14.413 

± 

0.048 

12.289 ± 0.394 

3.76 ± 0.30 

WISEA J140518.32+553421.3 

Y0.5 

18.765 ± 0.396 

14.097 

± 

0.037 

12.204 ± 0.263 

4.67 ± 0.40 

WISE J154151.65—225024.9 a 

Y1 

16.736 ± 0.165 b 

14.246 

± 

0.063 

>12.312 

2.49 ± 0.18 b 

WISEA J154214.00+223005.2 

T9.5 

18.846 ± 0.425 

15.043 

± 

0.061 

>13.014 

3.80 ± 0.43 

WISEA J 163940.84—684739.4 

YOpec 

17.266 ± 0.187 

13.544 

± 

0.059 

>11.755 

3.72 ± 0.20 

WISEA J173835.52+273258.8 

Y0 

17.710 ± 0.157 

14.497 

± 

0.043 

12.448 ± 0.399 

3.21 ± 0.16 

WISEA J205628.88+145953.6 

Y0 

16.480 ± 0.075 

13.839 

± 

0.037 

11.731 ± 0.249 

2.64 ± 0.08 

WISEA J220905.75+271143.6 

Y0: 

>18.831 

14.770 

± 

0.055 

12.455 ± 0.387 

>4.06 

WISEA J222055.34—362817.5 

Y0 

>18.772 

14.714 

± 

0.056 

>12.292 

>4.06 

Interlopers 

WISEA J013810.99+201657.5 


18.120 ± 0.273 

15.221 ± 0.087 

11.670 ± 0.199 

2.90 ± 0.29 

WISEA J065954.18—585559.6 


17.897 ± 0.133 

15.439 ± 0.059 

12.955 ± 0.371 

2.46 ± 0.15 


a WISE 1541 — 2250 is not in the A11WISE catalog. 
b Details of the likely erroneous W1 meas urement of WISE 1541 — 2250 
are discussed in |Kirkpatrick et al.| ( f2012| . 
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Table 2 

HST and Spitzer Aperture Photometry 


A11WISE Name 

F105W 

F125W 

F140W 

chi 

ch2 

chl—ch2 


(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 


WISEA J032504.52—504403.0 
WISEA J033515.07+431044.7 
WISEA J035000.31—565830.5 
WISEA J035934.07—540154.8 
WISEA J040443.50—642030.0 
WISEA J041022.75+150247.9 
WISEA J053516.87—750024.6 
WISEA J064723.24—623235.4 
WISEA J073444.03—715743.8 
WISEA J082507.37+280548.2 
WISEA J094306.00+360723.3 
WISEA J120604.25+840110.5 
WISEA J140518.32+553421.3 
WISE J154151.65—225024.9 
WISEA J154214.00+223005.2 
WISEA J163940.84—684739.4 
WISEA J173835.52+273258.8 
WISEA J205628.88+145953.6 
WISEA J220905.75+271143.6 
WISEA J221216.27—693121.6 
WISEA J222055.34—362817.5 
WISEA J235402.79+024014.1 


20.618 

± 

0.009 

19.598 

± 

0.002 

20.880 

± 

0.015 

20.092 

± 

0.003 

21.115 

± 

0.013 

20.276 

± 

0.003 

23.140 

± 

0.049 

22.801 

± 

0.051 

23.592 

± 

0.054 

23.453 

± 

0.050 

21.732 

± 

0.031 

20.964 

± 

0.004 

23.487 

± 

0.040 

23.197 

± 

0.030 

21.694 

± 

0.025 

21.062 

± 

0.005 

21.939 

± 

0.024 




22.204 

± 

0.044 

21.871 

± 

0.023 

21.268 

± 

0.020 

20.630 

± 

0.003 

21.337 

± 

0.026 

21.151 

± 

0.012 

23.842 

± 

0.057 




21.069 

± 

0.012 

20.347 

± 

0.003 

21.638 

± 

0.027 

20.997 

± 

0.005 




23.368 

± 

0.094 





17.746 

± 

0.086 




16.612 

± 

0.040 

22.321 

± 

0.047 

17.936 

± 

0.096 

21.806 

± 

0.039 

17.553 

± 

0.072 




17.633 

± 

0.082 

19.634 

± 

0.007 

16.636 

± 

0.042 

22.422 

± 

0.050 

17.753 

± 

0.084 




17.893 

± 

0.092 




17.649 

± 

0.077 




17.624 

± 

0.077 

20.037 

± 

0.007 

16.746 

± 

0.043 




17.339 

± 

0.061 

21.271 

± 

0.018 

16.876 

± 

0.046 




16.658 

± 

0.042 

20.240 

± 

0.010 

17.257 

± 

0.059 

19.883 

± 

0.008 

17.093 

± 

0.053 

19.524 

± 

0.007 

16.031 

± 

0.030 

23.167 

± 

0.149 

17.815 

± 

0.087 




17.364 

± 

0.063 




17.200 

± 

0.057 




18.105 

± 

0.109 


15.696 

± 

0.025 

2.050 

± 

0.090 

14.381 

± 

0.020 

2.231 

± 

0.045 

14.688 

± 

0.020 

3.248 

± 

0.098 

15.326 

± 

0.023 

2.227 

± 

0.075 

15.418 

± 

0.022 

2.216 

± 

0.085 

14.166 

± 

0.019 

2.470 

± 

0.046 

15.009 

± 

0.021 

2.744 

± 

0.087 

15.070 

± 

0.022 

2.823 

± 

0.094 

15.213 

± 

0.022 

2.436 

± 

0.080 

14.637 

± 

0.020 

2.987 

± 

0.079 

14.284 

± 

0.019 

2.461 

± 

0.047 

15.220 

± 

0.022 

2.119 

± 

0.065 

14.058 

± 

0.019 

2.818 

± 

0.050 

14.228 

± 

0.019 

2.430 

± 

0.046 

15.057 

± 

0.022 

2.200 

± 

0.062 

13.537 

± 

0.017 




14.473 

± 

0.019 

2.620 

± 

0.056 

13.923 

± 

0.018 

2.108 

± 

0.035 

14.739 

± 

0.020 

3.076 

± 

0.090 

14.973 

± 

0.021 

2.391 

± 

0.066 

14.736 

± 

0.021 

2.464 

± 

0.061 

15.013 

± 

0.022 

3.091 

± 

0.111 
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Table 3 

HST/WFC3 Spectroscopy Log 



A11WISE Name 

Date 

(UT) 

mode Total Int. Time 

(s) 

# of Exp. 

WISEA 

J032504.52—504403.0 

2013 July 30 

G102 

1612 

4 



2013 Aug 4 

G141 

1612 

4 

WISEA 

J033515.07+431044.7 

2013 July 12 

G102 

1812 

4 



2013 Aug 30 

G141 

1812 

4 

WISEA 

J035000.31—565830.5 a 

2011 Aug 13 

G141 

2212 

4 

WISEA 

J035934.07—540154.8 a 

2011 Aug 10 

G141 

2212 

4 

WISEA 

J040443.50—642030.0 

2013 April 9 

G102 

1812 

4 



2013 April 9 

G141 

1812 

4 

WISEA 

J041022.75+150247.9 

2012 Sep 1 

G141 

2012 

4 

WISEA 

J053516.87—750024.6 a 

2011 Sep 27 

G141 

2212 

4 



2012 Sep 17 

G141 

2212 

4 



2013 Sep 26 

G102 

7618 

6 



2013 Sep 27 

G141 

7615 

6 



2013 Sep 27 

G102 

7618 

6 



2013 Dec 4 

G102 

7618 

6 

WISEA 

J064723.24—623235.4 b 

2013 May 13 

G141 

7218 

6 



2013 May 13 

G102 

7036 

12 



2013 May 15 

G102 

7036 

12 



2013 Nov 14 

G102 

7036 

12 



2013 Dec 28 

G102 

7036 

12 



2013 Dec 29 

G141 

7218 

6 



2013 Dec 30 

G102 

7036 

12 

WISEA 

J073444.03—715743.8 

2013 May 18 

G102 

1812 

4 



2013 May 20 

G141 

1812 

4 

WISEA 

J082507.37+280548.2 

2014 Jan 16 

G102 

7018 

6 



2014 Jan 17 

G141 

7218 

3 



2014 Jan 18 

G102 

7018 

6 



2014 Jan 19 

G102 

7018 

6 

WISEA 

J094306.00+360723.3 

2013 Feb 20 

G141 

2012 

4 

WISEA 

J120604.25+840110.5 

2013 July 15 

G102 

1812 

4 



2013 July 15 

G141 

1812 

4 

WISEA 

J140518.32+553421.3° 

2011 Mar 14 

G141 

2212 

4 



2013 Apr 18 

G102 

1812 

4 

WISE J154151.65—225024.9 

2013 May 9 

G102 

4612 

4 



2013 May 9 

G141 

1812 

4 

WISEA 

J154214.00+223005.2 

2012 Mar 4 

G141 

2012 

4 



2013 June 5 

G102 

1812 

4 



2013 June 7 

G141 

1812 

4 

WISEA 

J 163940.84—684739.4 

2013 Oct 26 

G102 

7518 

6 



2013 Oct 27 

G102 

7518 

6 



2013 Oct 29 

G141 

10024 

4 

WISEA 

J 173835.52+273258.8° 

2011 May 12 

G141 

2012 

4 

WISEA 

J205628.88+145953.6 a 

2011 Sept 4 

G141 

2012 

4 

WISEA 

J220905.75+271143.6 d 

2012 Sept 15 

G141 

2012 

4 



2013 Apr 28 

G102 

7218 

6 



2013 June 6 

G102 

7218 

6 



2013 Sept 20 

G102 

7218 

6 

WISEA 

J221216.27-693121.6 

2013 Sept 7 

G102 

4612 

4 



2013 Sept 11 

G141 

1812 

4 

WISEA 

J222055.34—362817.5 

2013 June 8 

G102 

1812 

4 



2013 June 20 

G141 

4412 

4 

WISEA 

J235402.79+024014.1 

2013 Sept 22 

G141 

3224 

8 


a Published previously in 
b Published previously in 
c Published previously in 
d Published previously in 


KirkjDatricl^^+ (2012 


Kirkpatrick et al. (2013 


fcT^h et al. 


( 20111 ). 


Cushing et al 


2014aJ. 
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Table 4 

Distance Estimates 


A11WISE Name 

Spec. 

Dist a 

Dist b 


Type 

(pc) 

(pc) 

WISEA J032504.52—504403.0 

T8 

36.4+2.4 

36.0+2.4 

WISEA J040443.50—642030.0 

T9 

24.8+0.7 

24.5+0.7 

WISEA J082507.37+280548.2 

Y0.5 

10.0+0.3 

10.9+0.3 

WISEA J120604.25+840110.5 

Y0 

14.5+0.4 

15.0+0.4 

WISEA J221216.27-693121.6 

T9 

16.8+0.5 

16.6+0.5 

WISEA J235402.79+024014.1 

Y1 

10.1+0.4 

11.9+0.5 


a Using W2 relation of Kirkpatrick et al. 1|2Q12| ). 
b Using W2 relation of Dupu^^^Liuf|2Q12| ) . 


Table 5 

Synthetic Photometry 


A11WISE Name 

Spec. 

Ymko 

J MKO 

H MKO 

F105W 

F125W 

F140W 


Type 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 


WISEA 

J032504.52—504403.0 

T8 

19.980 

± 

0.027 

18.935 

± 

0.024 

19.423 

± 

0.027 

20.601 

± 

0.028 

19.547 

± 

0.023 

19.223 

± 

0.022 

WISEA 

J033515.07+431044.7 

T9 

20.166 

± 

0.029 

19.467 

± 

0.023 

19.938 

± 

0.031 

20.939 

± 

0.033 

20.137 

± 

0.025 

19.785 

± 

0.023 

WISEA 

J035000.31—565830.5 a 

Y1 




22.178 

± 

0.073 

22.263 

± 

0.135 




22.951 

± 

0.114 

22.431 

± 

0.073 

WISEA 

J035934.07—540154.8 a 

Y0 




21.566 

± 

0.046 

22.028 

± 

0.112 




22.258 

± 

0.062 

21.789 

± 

0.045 

WISEA 

J040443.50—642030.0 

T9 

20.328 

± 

0.032 

19.647 

± 

0.025 

19.970 

± 

0.033 

21.063 

± 

0.037 

20.293 

± 

0.028 

19.893 

± 

0.024 

WISEA 

J041022.75+150247.9 a 

Y0 




19.325 

± 

0.024 

19.897 

± 

0.038 




19.997 

± 

0.025 

19.643 

± 

0.024 

WISEA 

J053516.87—750024.6 b 

>Y1 

22.701 

± 

0.070 

22.132 

± 

0.071 




23.581 

± 

0.138 

22.876 

± 

0.102 




WISEA 

J064723.24—623235.4 

Y1 

22.870 

± 

0.076 

22.854 

± 

0.066 

23.306 

± 

0.166 

23.833 

± 

0.117 

23.683 

± 

0.098 

23.204 

± 

0.066 

WISEA 

J073444.03—715743.8 

Y0 

20.870 

± 

0.041 

20.354 

± 

0.029 

21.069 

± 

0.071 

21.675 

± 

0.051 

21.045 

± 

0.035 

20.726 

± 

0.030 

WISEA 

J082507.37+280548.2 

Y0.5 

22.566 

± 

0.053 

22.401 

± 

0.050 

22.965 

± 

0.139 

23.409 

± 

0.073 

23.015 

± 

0.062 

22.731 

± 

0.051 

WISEA 

J094306.00+360723.3 a 

T9.5 




19.766 

± 

0.025 

20.315 

± 

0.038 




20.444 

± 

0.027 

20.092 

± 

0.025 

WISEA 

J120604.25+840110.5 

Y0 

20.875 

± 

0.036 

20.472 

± 

0.030 

21.061 

± 

0.062 

21.819 

± 

0.050 

21.171 

± 

0.036 

20.798 

± 

0.029 

WISEA 

J140518.32+553421.3 

Y0.5 

21.333 

± 

0.057 

21.061 

± 

0.035 

21.501 

± 

0.073 

22.193 

± 

0.076 

21.730 

± 

0.044 

21.375 

± 

0.035 

WISE J154151.65—225024.9 

Y1 

21.671 

± 

0.037 

21.631 

± 

0.064 

22.085 

± 

0.170 

22.724 

± 

0.083 

22.443 

± 

0.097 

22.138 

± 

0.078 

WISEA 

J154214.00+223005.2 

T9.5 

20.461 

± 

0.028 

19.937 

± 

0.026 

20.520 

± 

0.045 

21.273 

± 

0.034 

20.586 

± 

0.027 

20.235 

± 

0.025 

WISEA 

J 163940.84—684739.4 

YOpec 

20.833 

± 

0.023 

20.626 

± 

0.023 

20.746 

± 

0.029 

21.674 

± 

0.025 

21.252 

± 

0.024 

20.824 

± 

0.024 

WISEA 

J173835.52+273258.8 a 

Y0 




19.546 

± 

0.023 

20.246 

± 

0.031 




20.223 

± 

0.023 

19.923 

± 

0.023 

WISEA 

J205628.88+145953.6 a 

Y0 




19.129 

± 

0.022 

19.643 

± 

0.026 




19.811 

± 

0.023 

19.479 

± 

0.022 

WISEA 

J220905.75+271143.6 

Y0: 

22.954 

± 

0.071 

22.859 

± 

0.128 

22.389 

± 

0.152 

23.508 

± 

0.110 

23.355 

± 

0.156 

22.877 

± 

0.105 

WISEA 

J221216.27-693121.6 

T9 

20.282 

± 

0.023 

19.737 

± 

0.024 

20.225 

± 

0.036 

21.043 

± 

0.027 

20.378 

± 

0.026 

20.047 

± 

0.024 

WISEA 

J222055.34—362817.5 

Y0 

20.899 

± 

0.034 

20.447 

± 

0.025 

20.858 

± 

0.035 

21.783 

± 

0.042 

21.131 

± 

0.027 

20.749 

± 

0.024 

WISEA 

J235402.79+024014. l a 

Y1 




23.068 

± 

0.199 

22.882 

± 

0.300 




24.124 

± 

0.393 

23.292 

± 

0.179 


a Y-band and F105W photometry were not synthesized for these dwarfs 
because they were not observed with the G102 grism. 
b H-band and F140W photometry were not synthesized for WISE 
0535—7500 because its spectrum is contaminated in this wavelength 
range. 
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Table 6 

Models Used for Spectral Fitting® 


T eff 

(K) 

log 9 
(cm s -1 ) 

fsed 

Cloud Type 
(%) 

Reference 

400-700(50), 700-1000(100) 

3.0 

nc 


1 

400-1000(50) 

3.5 

nc 


1 

300-1000(50) 

3.75 

nc 


1 

300-1000(50) 

4.0 

nc 


1 

400-1000(50) 

4.25 

nc 


1 

300-1000(50) 

4.5 

nc 


1 

400-1000(50) 

4.75 

nc 


1 

400-1000(50) 

5.0 

nc 


1 

400-1000(100) 

5.5 

nc 


1 

200-400(25) 

3.0 

5 

50 b 

2 

200-400(25), 450 

3.5 

5 

50 b 

2 

200-400(25), 450 

4.0 

3 

50 b 

2 

200-400(25), 450 

4.0 

5 

50 b 

2 

200-400(25), 450 

4.0 

7 

50 b 

2 

200-400(25), 450 

4.5 

5 

50 b 

2 

200-400(25), 400-700(50), 700-1000(100) 

5.0 

5 

50 b 

2 

400, 450, 550, 600-1000(100) 

4.0 

2 

100 

3 

400-600(50), 600-1000(100) 

4.0 

3 

100 

3 

400-600(50), 600-1000(100) 

4.0 

4 

100 

3 

400-600(50), 600-1000(100) 

4.0 

5 

100 

3 

400-600(50), 600-1000(100) 

4.5 

2 

100 

3 

400-600(50), 600-1000(100) 

4.5 

3 

100 

3 

400-600(50), 600-1000(100) 

4.5 

4 

100 

3 

400-600(50), 600-1000(100) 

4.5 

5 

100 

3 

460, 500-600(50), 600-1000(100) 

5.0 

2 

100 

3 

400-600(50), 460, 600-1000(100) 

5.0 

3 

100 

3 

400, 460, 500-600(50), 600-1000(100) 

5.0 

4 

100 

3 

400, 460, 500-600(50), 600-1000(100) 

5.0 

5 

100 

3 

400, 500-600(50), 600- 1000(100) 

5.5 

3 

100 

3 


References: (1) Saumon et al.|2012| (2)|Morley et al.|2014| (3)|Morley| 
let al. 120121 

a When a range is given, the value given in parentheses is the spacing 
for that interval. 
b With water clouds. 
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Table 7 

Model-derived Physical Characteristics 3, 


A11WISE Name 

Spec. 

Type 

T eff 

(K) 

log 9 

(cm s -1 ) 

fsed 

% Cloudy 
(%) ' 

(R/d) 
O&jup/ P c ) 

X 2 /d.o.f 

Age b 

(Gyr) 

WISEA J032504.52—504403.0 

T8 

550-600 

4.0 

3-4 

100 

2.468xl0 -2 

5213/211 

0.08-0.3 

WISEA J033515.07+431044.7 

T9 

500-550 

4.0—4.5 

3-5 

100 

1.973xl0 -2 

5276/211 

0 .2-1.5 

WISEA J035000.31—565830.5 

Y1 

300-350 

5.0 

... ,5 

0,50 

4.704xl0 -2 

2873/128 

U 

WISEA J035934.07—540154.8 

Y0 

400 

4.0—5.0 

2,3,5 

100 

3.641X10 -2 

2494/128 

1.5-8 

WISEA J040443.50—642030.0 

T9 

550-600 

4.0 

3-4 

100 

1.861X10 -2 

4458/211 

0.08-0.3 

WISEA J041022.75+150247.9 

Y0 

400 

4.0—4.5 

2,3,4 

100 

1.007xl0 _1 

5172/128 

1.5-8 

WISEA J053516.87—750024.6 

>Y1 

450-500 

5.0 

... ,5 

0,50 

8.397xl0 -3 

2744/179 

4—U 

WISEA J064723.24—623235.4 

Y1 

350-400 

5.0 

... ,5 

0,50 

2.609xl0 -2 

2621/211 

8 —U 

WISEA J073444.03—715743.8 

Y0 

450 

4.0—4.5 

3-5 

100 

2.410xl0 -2 

3109/211 

0.4-2 

WISEA J082507.37+280548.2 

Y0.5 

400 

4.5—5.0 

... ,4,5 

0,100 

1.808xl0 -2 

3313/211 

3—U 

WISEA J094306.00+360723.3 

T9.5 

450-500 

4.0—4.5 

2,3,5 

100 

4.167xl0 -2 

4506/128 

0.3-2 

WISEA J120604.25+840110.5 

Y0 

400-450 

4.0—4.5 

3,5 

100 

4.040xl0 -2 

3475/211 

0.4-3 

WISEA J140518.32+553421.3 

Y0.5 

350-400 

5.0—5.5 

... ,5 

0,100 

3.218xl0 -2 

4117/211 

U 

WISE J154151.65-225024.9 

Y1 

400 

4.0—4.5 


0 

2.072xl0 -2 

4102/211 

0.6-3 

WISEA J154214.00+223005.2 

T9.5 

450-500 

4.0—4.5 

3-5 

100 

2.127xl0 -2 

4051/211 

0.3-2 

WISEA J 163940.84—684739.4 

YOpec 

400 

5.0 

... ,5 

0,100 

3.186xl0 -2 

10859/211 

8 —U 

WISEA J 173835.52+273258.8 

Y0 

400 

4.0—4.5 

2,4,5 

100 

6.445xl0 -2 

6598/128 

0.6-3 

WISEA J205628.88+145953.6 

Y0 

400-450 

4.0—4.5 

2,4 

100 

5.541xl0 -2 

10770/128 

0.4-3 

WISEA J220905.75+271143.6 

Y0: 

500-550 

4.0—4.5 

4,5 

100 

6.682xl0 -3 

3106/211 

0 .2-1.5 

WISEA J221216.27—693121.6 

T9 

500-600 

4.0 

3,5 

100 

1.541xl0 -2 

5389/211 

0.08-0.4 

WISEA J222055.34—362817.5 

Y0 

400-450 

4.0—5.0 

2,3,5 

100 

3.204xl0 -2 

4574/211 

1-6 

WISEA J235402.79+024014.1 

Y1 

300-400 

4.0—5.0 

2,5 

50,100 

5.961X10 -2 

10663/128 

1.5—U 


a Model-derived parameters are generally unreliable for reasons ex¬ 
plained in Section 5.4 

b Ages are derived for the T e ff and log g ranges in columns 3 and 
4 using the Saumon & Marley (2008) cloudless evolutionary models. 
Some combinations of T e f / and log g derived from the atmospheric 
fits are unphysical in that brown dwarfs cannot evolve to have such 
values, so the corresponding ages are denoted as ‘U’ for unphysical. 











